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基于 线粒体 COI 基因 序列 的 大 豆 食 心 虫 
中 国 东 北 地 理 种 群 遗 传 多 样 性 分 析 
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摘要 :【 目 的 ] 大 豆 食心虫 Keguminizora glycinivorella (Matsumura) 是 一 种 危害 大 豆 的 主要 害虫 ,在 中 国 北 方 地 区 危 
害 较 重 。 本 研究 旨 在 探讨 大 豆 食心虫 在 中 国 东北 不 同 地 理 种 群 间 的 遗传 变异 。[ 方 法 ] 测 定 了 10 个 不 同 地 理 种 群 
153 个 个 体 的 线粒体 细胞 色素 氧化 酶 亚 基 工 (mtCOI) 基因 的 657 bp 序列 ,利用 DnaSP 5. 0 和 Arlequin 3. 5. 1. 2 
等 软件 对 大 豆 食心虫 种 群 间 的 遗传 多 样 性 .基因 流水 平和 分 子 变异 进行 分 析 。【 结 果 ]】 结 果 表 明 :10 个 地 理 种 群 间 
的 COI 基因 共有 36 个 变异 位 点 和 17 个 单 倍 型 ,其 中 1 个 单 倍 型 为 10 个 种 群 所 共享 。 总 种 群 的 单 倍 型 多 样 性 指数 
Hd 为 0.456, 各 地 理 种 群 单 倍 型 多 样 度 范围 在 0 ~0.634 之 间 。 总 群体 的 固定 系数 Fst 为 0. 12545 , 遗传 分 化 系数 
Gst 为 0.06326 ,总 基因 流 Nm 为 3.49, 且 各 种 群 间 的 基因 流 均 大 于 1, 种 群 间 基因 交流 的 水 平 较 高 。【 结论 ] 大 豆 食 
心 忠 种 群 内 遗传 多 样 性 水 平 处 于 中 低 等 水 平 。 总 群体 和 各 种 群 的 Tajima' s D 检验 结果 皆 不 显著 ,说 明 中 国 东 北 地 
区 大 豆 食 心 虫 在 较 近 的 历史 时 期 内 没有 出 现 种 群 扩张 现象 。 AMOVA 分 子 变异 分 析 结 果 表 明 , 大 豆 食心虫 的 遗传 
分 化 主要 来 自 种 群 内 部 ,而 种 群 间 未 发 生 明 显 的 遗传 分 化 。 各 地 理 种 群 的 单 倍 型 在 系统 发 育 树 上 和 中 介 网 络 图 上 
散布 在 不 同 的 分 布 群 中 ,缺乏 明显 的 地 理 分 布 格局 。 各 种 群 的 遗传 距离 与 地 理 距 离 之 间 没 有 显著 线性 相关 性 ,种 
群 间 的 基因 交流 并 未 受到 地 理 距 离 的 影响 。 

关键 词 : 大 豆 食心虫 ; 线粒体 COI 基因 ; 单 倍 型 ; 地 理 种 群 ; 遗传 多 样 性 ; 基因 流 
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Analysis of the genetic diversity in geographic populations of Leguminivora 
glycinivorella ( Lepidoptera: Olethreutidae) from northeastern China based 


on mitochondrial DNA COI gene sequences 
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Agronomy, Northeast Agricultural University, Harbin 150030, China; 2. Institute of Plant Protection, 
























































































































































































































































Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China) 

Abstract: [ Aim] The soybean podborer, Leguminivora glycinivorella ( Matsumura) , is a major insect 
pest of soybean Glycina max, and it damages G. max seriously in northeastern China. This study aims to 
investigate the genetic differentiation among geographic populations of L. glycinivorella in northeastern 
China. [Methods] A 657-bp segment of the mitochondrial cytochrome c oxidase subunit I ( COI) gene of 
153 individuals from 10 geographic populations in northeastern China was sequenced and analyzed. The 
genetic diversity, gene flow and molecular variance were analyzed using DnaSP 5. 0 and Arlequin 
3. 5. 1. 2. [Results] The results showed that 36 variable sites were detected, and 17 haplotypes were 
defined based on the COI sequences of 10 geographic populations of L. glycinivorella, including one 
haplotype shared by all the 10 populations. The total haplotype diversity of all populations was 0. 456, 
and the range of haplotype diversity of different populations was 0 — 0. 634. The total Fst, Gst and Nm 
were 0. 12545 , 0. 06326 and 3.49 , respectively, and all pairwise gene flow Nm was greater than 1 in the 
10 populations, suggesting that extensive gene flow occurred among the ten populations of L. 
glycinivorella. [ Conclusion] The L. glycinivorella populations are characterized by medium/low genetic 


diversity. The Tajima' s D test indicated that there might not have been a recent population expansion. 
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Molecular variance analysis ( AMOVA) demonstrated that the observed genetic differentiation occurs 


primarily within populations, and there is no obvious genetic differentiation among populations. The 


neighbor-joining phylogenetic tree and haplotype network showed that the haplotypes are distributed in 


different clades and no obvious geographical structure has been formed. The genetic distance is not 


significantly correlated with the geographical distance between populations, and the gene flow is not 


affected by geographical distance. 


Key words: Leguminivora glycinivorella; mtDNA COI gene; haplotype; geographic population; genetic 


diversity; gene flow 


AK Ea B Leguminivora glycinivorella ( Mats. ) 
XU CK xi RE SE SR. UN EL Hh. 豆荚 虫 , 属 鳞 翅 H 
( Lepidoptera) ,小 卷 峨 科 ( Olethreutidae) 。 在 中 国 主 
要 分 布 于 东北 ,华北 ,西北 及 湖北 、 江 苏 浙江 安徽 、 
山东 等 地 ,以 东北 三 省 受害 最 为 严重 , 日 本 、 朝 鲜 及 
远东 沿海 边区 也 有 分 布 。 其 食性 相对 较 单 一 ,以 幼 
虫 钻 性 咬 食 豆 粒 为 害 ( 杨 啸 等 ,2012 ) , 主要 为 害 大 
豆 Glycina max , 野生 大 豆 Glycine ussuriensis 和 和 苦 参 
Sopha flavescens, Æ F H, 虫 食 率 通常 为 10% ~ 
30% ,严重 时 可 达 50% 以 上 , 因 地 区 、 年 度 气 候 变 
化 大 豆 品种 的 不 同 而 差别 很 大 ,严重 影响 大 豆 产量 
和 品质 ( 匡 炳 芝 ,2001 ) 。 

目前 尚 没有 从 遗传 分 化 的 分 子 特 性 对 大 豆 食 
心 虫 地 理 种 群 进行 遗传 变异 分 析 的 相关 报道 。 随 
着 分 子 生物 学 的 发 展 和 应 用 ,分 子 标记 越 来 越 多 
地 被 应 用 于 昆虫 遗传 研究 中 去 ,在 众多 的 分 子 标 
记 技 术 中 ,线粒体 DNA 具有 提取 技术 简单 .母系 
遗传 .进化 速率 较 核 基因 人 快 等 特点 (Harrison， 
1989) ,被 广泛 应 用 于 昆虫 系统 进化 及 种 群 的 遗传 
变异 中 (Arvise et al., 1987; Barrio et al., 1992; 
Miya et al., 2003) 。 例 如 ,COL,COI 和 168 rDNA 
等 经 常用 于 分 析 亲 缘 关 系 密 切 的 种 、 亚 种 及 地 理 
种 群 之 间 的 系统 关系 (Xu Q et al., 2009; FE, 
2010; 29 B Æ 55,2011; Xu DD et al., 2012; $^ 
56,2013) 。 本 研究 针对 大 豆 食 心 虫 因 迁 飞 能 力 弱 
而 导致 长 期 地 理 隔离 的 特点 ,利用 线粒体 COT 基 
因 作 为 分 子 标记 ,测定 了 东北 地 区 大 豆 食 心 虫 10 
个 地 理 种 群 线粒体 COI 基因 序列 ,以 胡桃 小 卷 蛾 
Cydia caryana 的 线粒体 COI 基因 片段 (658 bp) 
( GenBank 登录 号 :GU088344. 1) 为 对 照 , 从 种 群 
内 及 种 群 间 的 基因 遗传 分 化 程度 和 基因 交流 水 平 
对 大 豆 食 心 虫 的 地 理 种 群 遗 传 多 样 性 进行 了 分 
析 ,为 探讨 大 豆 食 心 虫 不 同 地理 种 群 间 的 系统 发 
育 及 进化 关系 提供 遗传 学 信息 ,进而 为 该 虫 的 区 
域 性 控制 和 防治 策略 提供 依据 。 
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1 材料 和 方法 


1.1 供 试 虫 源 

本 研究 以 东北 三 省 大 豆 食 心 虫 危 害 较 重 的 10 
个 地 理 种 群 的 老 熟 幼虫 作为 供 试 虫 源 。 寄 主 均 为 当 
地 种 植 的 栽培 大 豆 , 采 用 随机 取样 的 方法 进行 采集 ， 
在 大 豆 晚 熟 期 拨 开 表面 有 大 豆 食 心 由 外 性 痕迹 的 豆 
美 ,用 毛笔 挑 出 老 熟 幼虫 ,将 其 放 入 成 有 无 水 乙醇 的 
离心 管 中 ,标明 采集 样品 信息 , - 207€ 的 冰箱 中 保存 
备用 。 具 体 样 品 采集 的 地 点 、 经 纬度 海拔 等 信息 详 
见 表 1。 
1.2 基因 组 DNA 的 提取 

大 豆 食心虫 单 头 老 熟 幼 虫 经 双 芝 水 漂洗 、 吸 水 
纸 吸 干 DNA 裂解 液 (100 mmol/L Tris-HCl pH 8.0, 
25 mmol/L EDTA ,500 mmol/L NaCl ,1% SDS) 裂解 ， 
之 后 用 酚 氯 仿 抽 提 法 提取 虫 体 基 因 组 DNA 
(Sambrook and Russell ，2001 ) 。 提 取 的 DNA 样品 经 
琼脂 糖 凝 胶 电泳 检测 后 保存 于 -20%C 冰 箱 中 备用 。 
1.5 目的 片段 PCR 扩 增 及 电泳 

本 研究 线粒体 COL PCR 所 用 的 引物 参照 
Folmer && (1994) 的 通用 引物 LCO1490: 5'-GGT 
CAACAAATCATAAAGATATTGG-3' 和 HCO2198;5'- 
TAAACTTCAGGGTGACCAAAAAATCA-3', PCR X 
应 体系 为 50 uL, HEHP 10 x PCR buffer 5 uL, dNTPs 4 
uL,DNA 模板 3 pL,0.01 mmol/L 的 上 游 和 下 游 引物 
各 1 uL, Taq DNA 聚合 酶 ( 5 U/uL) 0.5 mL。 扩 增 条 
件 为 94% 预 变 性 5 min; 然 后 94% 变性 45 s,53 CB c 
50 s,72% 延伸 45 s, 共 35 个 循环 ;72% 后 延伸 10 
min,4% 保存 。PCR 产物 采用 196 琼脂 糖 凝 胶 检 测 ， 
然后 在 凝 胶 成 像 分 析 系 统 中 观察 并 照相 ,记录 结果 。 
1.4 PCR 产物 的 纯化 回收 及 测序 

按照 DNA 凝 胶 回 收 试剂 盒 说 明 书 回收 各 样品 
PCR 产物 。 回 收 产物 委托 北京 华 大 基因 研究 中 心 
进行 DNA 测序 。 
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表 1 栽培 大 豆 为 寄主 的 大 豆 食心虫 不 同 地 理 种 群 样 本 信息 


Table 1 Specimen data of difference geographic populations of Leguminivora glycinivorella on cultivated soybean 





























种 群 代码 样本 数 采集 地 点 AEE) 纬度 CN) 海拔 (my) 采集 时 间 
Population code Number of specimens Collecting locality Longitude Latitude Elevation Collecting date 
i^» 
SY 17 E 123734! 41°49" 51 2013. 10 
Shenyang, Liaoning 
Ba dry 哈 尔 
QQHR 15 黑龙 江 齐 齐 b 123555' 47?22' 149 2013.9 
Qiqihar, Heilongjiang 
去 林 公 主 岭 
GZL 15 z 林 公 id 124?49' 43?30' 212 2013.10 
Gongzhuling, Jilin 
黑龙 江 大 庆 
DQ 13 黑龙 江 大 有 a 125?*10" 46?33' 149 2013.9 
Daqing, Heilongjiang 
吉林 德 惠 
DH 10 EE 125*57' 44?4]" 174 2013.10 
Dehui, Jilin 
黑龙 江 哈 尔 滨 
HRB 18 黑龙 江 TS 126?45' 45?43' 118 2013.9 
Harbin, Heilongjiang 
Ba yT 
SH 12 ] 黑龙 江 绥化 126*56' 46?37' 175 2013.9 
Suihua, Heilongjiang 
58 JE TE SER 
HH 24 黑龙 江 黑 河 127*28' 50?15' 122 2013.9 
Heihe, Heilongjiang 
PIT y 
MDJ 12 黑龙 江 牡丹江 a 129?30' 44?25' 230 2013.9 
Mudanjiang, Heilongjiang 
EVA 
JMS 17 黑龙 江 佳木斯 130?28' 46^49' 82 2013.9 


Jiamusi, Heilongjiang 


1.5 数据 统计 与 序列 分 析 方法 
应 用 分 子 进化 遗传 分 析 软 件 MEGA. 5. 0 
,2011) 对 DNA 序列 进行 序列 比 对 及 
同 源 性 分 析 ， de USD 保守 位 点 、 
变异 位 点 、 简 约 信息 位 点 、 单 一 突变 位 点 、 转 换 碱 基 
对 、 赤 换 碱 基 对 数 、 转 换 和 匡 换 的 比值 以 及 单 倍 型 间 
的 遗传 距离 。 利 用 DnaSP 5. O 软件 ( Librado and 
Rozas, 2009 ) 统计 COIL 基因 序列 单 倍 型 数目 和 出 现 
的 频率 、 单 倍 型 多 样 度 ( Hd) AEREE CPi) \ 核 
苷 酸 平均 差异 数 (K) 、 种 群 间 遗 传 分 化 程度 参数 Fst 
和 Gst ,并 进行 中 性 检验 (Tajima，1989 ) 。 根 据 线 粒 
体 基 因 特 有 公式 Fst 21/(1 +2Nm) (Takahata and 
Palumbi, 1985; Goldberg and Ruvolo, 1997) ,推算 基 
因 流 Nm 值 。 根 据 经 纬度 计算 出 各 采集 地 之 间 的 地 
理 距 离 。 采 用 Mega 5. 0 对 序列 进行 系统 发 育 分 析 ， 
根据 Kimura 2-Parameter 选用 邻接 法 (NJ) 
( Felsenstein, 1985) 构建 大 豆 食心虫 不 同 地 理 种 群 
的 分 子 系统 进化 树 ,进化 树 分 支 的 置信 度 采 用 自 展 
法 (Bootstrap ) 重复 检测 1 000 次 。 应 用 Arlequin 
3.5.1.2 软 件 ( Excoffier et al., 2005) ,根据 Pairwise 
difference 模型 ,计算 种 群 间 的 分 化 指数 (Fst) 和 分 
子 方 差分 析 以 评估 种 群 之 间 、 种 群 内 的 遗传 分 化 。 
采用 Network 4. 6 程序 ( Bandelt et al., 1999) 构建 种 

群 内 所 有 单 倍 型 的 中 介 网 络 图 。 


(Tamura et al. 





























2 结果 与 分 析 


2.1 大 豆 食 心 虫 COI 基因 的 碱 基 组 成 及 序列 分 析 

以 基因 组 DNA 为 模板 进行 PCR 扩 增 ,从 所 有 
测序 的 样品 中 均 得 到 一 条 特异 的 PCR 产物 ,最 后 得 
到 657 bp 的 COI 基因 序列 ( 表 2)。 测 序 结果 在 
NCBI 中 用 BLAST 进行 相似 性 检索 ,与 胡桃 小 卷 峨 
Cydia caryana COI 基因 序列 (GenBank 登录 号 : 
GU088344. 1) 的 核 昔 酸 一 致 性 为 93% ,为 所 需 目的 
片段 。 

本 研究 共 获 得 153 条 序列 , 均 没 有 碱 基 的 缺失 
HIA. Hr IH AR Ceonserved sites) 621 个 , 变 
异 位 点 (variable sites) 36 个 ( 表 2) ,变异 百分率 为 
5.4896 ,突变 总 数 (total number of mutation ) 为 37。 
其 中 简约 信息 位 点 (parsimony informative sites ) 24 
个 ,单一 变异 位 点 (singleton variable sites ) 12 个 。 i 
测序 列 中 COL 基因 碱 基 组 成 为 :A 30. 396, 
37.896 ,C 17.0% 和 G14.9% ,A+T 含 量 为 68.12% , 
G+C REN 31.9% ,表现 出 明显 的 碱 基 AT f fg 
性 ,符合 昆虫 线粒体 基因 序列 碱 基 组 成 特点 (Jermiin 
and Crozier, 1994) 。 所 有 的 核 苷 酸 替 换 分 析 表 明 ， 
转换 (transition ) 占 88. 89% ,其 中 A 5 G 之 间 的 转 
换 占 22.22% ,C 与 T 之 间 的 转换 占 66. 67% ; Bj d 
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(transversion ) 发 生 频 率 较 转 换 低 为 11. 1196 ;总 体 — BOATS TIC IRVECIT RRR ARRAK T BR 
变异 位 点 中 转换 与 颠 换 的 比值 为 8:1 ,与 亲缘 关系 ” 换 的 变异 特点 (Simon et al., 1994) 一 致 。 





R2 大 豆 食心虫 不 同 单 倍 型 的 COI 基因 序列 变异 
Table 2 Sequence variations of different haplotypes of COI gene in Leguminivora glycinivorella 
变异 位 点 Variable sites 
2 2 2 
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变异 位 点 Variable sites 
5 5 5 





单 倍 型 
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H13 . . . . í : l : . : i . A : C 
H14 
H15 
H16 : : à f ; ; l ; ; T 
H17 . : . . : : f G : f f : : 
序列 全 长 为 657 bp ,对 齐 后 有 36 个 变异 位 点 。HIl ~ HIT 为 大 豆 食心虫 的 07 个 单 倍 型 ,没有 碱 基 的 缺失 和 插入 。 点 号 表示 相同 的 碱 基 。Thirty 


six variable sites were defined from the alignments with a total of 657 bp in length. H1 — H17 refer to the 17 haplotypes of L. glycinivorella ; no deletion/ 
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2.2 大豆 食 心 虫 不 同 地 理 种 群 COI 基因 的 单 倍 
型 、 核 苷 酸 多 样 度 和 中 性 检验 分 析 

上 述 153 条 COL 基因 序列 检测 出 17 种 单 倍 型 ， 
分 别 命名 为 H1 ~ H17( GenBank 登录 号 :KJ528252 ~ 
KJ528258 ;KJ540168 ~ KJ540177 ) 。 对 各 种 群 单 倍 型 
数 进行 多 样 性 分 析 和 Tajima' s D 41r C 3) RHH , f. 
倍 型 Hl 出 现 频率 最 高 为 73.2% (112/153) , H% 10 
个 地 理 种 群 所 共享 ,而 牡丹 江 ( MDJ) 的 全 部 个 体 均 属 
于 该 单 倍 型 。 其 次 是 单 倍 型 H2 出 现 的 频率 较 高 ， 
约 占 8. 5% (13/153) ,为 4 个 地 理 种 群 所 共有 。 
H12 有 56 个 个 体 ,H4 有 4 头 个 体 ,H5 有 2 头 个 体 ， 
H3 有 3 头 个 体 ,H6 有 3 头 个 体 ,剩余 其 他 单 倍 型 仅 
有 一 个 个 体 ,属于 稀有 单 倍 型 。 

















根据 DnaSP5. 0 软件 分 析 结果 表明 ,总 种 群 单 
倍 型 多 样 度 (haplotype diversity, Hd) 为 0.456 , zd 
酸 多 样 度 (nucleotide diversity, Pi) 为 0. 00676 
(Smith et al., 2006) , 4% P IRE 35] 25 EC ( average 
number of nucleotide differences, K) X 4. 440, 不 同 
地 理 种 群 的 单 倍 型 多 样 度 在 0 ~0. 634 zt [8] , ECHEBR 
多 样 度 在 0 ~0.01435 之 间 ( 表 3 ) ,其 中 牡丹 江 种 群 
各 项 参数 均 为 0。 从 Tajima’ s D 检测 结果 可 以 得 出 
(3€ 3) ,10 个 种 群 的 中 性 检验 结果 均 不 显著 , 且 总 
种 群 的 Tajima’ s D 值 为 - 0. 97226, 结果 不 显著 
(P >0. 10) ( Harpending et al., 1998) ,说 明 中 国 东 
北三 省 地 区 大 豆 食心虫 总 体 上 在 较 近 的 历史 时 期 未 
经 历 明 显 的 种 群 扩 张 。 





























X3 大 豆 食心虫 不 同 地 理 种 群 COI 基因 单 倍 型 多 样 度 、 核 昔 酸 多 样 性 分 析 及 Tajima' s D 中 性 检验 
Table 3 COI haplotype diversity, nucleotide diversity and Tajima’ s D test in different 


geographic populations of Leguminivora glycinivorella 





n 





PARNA PENET 








核 音 酸 多 样 度 


核 董 酸 平均 差异 数 








羊 代码 性 检验 性 检测 显著 性 
P ul d Haplotype distribution Haplotype diversity Nucleotide diversity Average number of Eid "à 人 性 
opulation code (h) (Hd) (Pi) nucleotide differences (K) ajima' s Statistical significance 
DQ H1(12) H7(1) 0.154 0.00047 0.308 — 1.46802 P »0.10 
DH H1(8) H9(1) HIO(1) 0.378 0.00237 1.556 - 1.57285 0.10 >P 20.05 
GZL H1(11) H2(2) H5(2) 0.457 0.01377 9.048 1. 39133 P »0.10 
HEB SCORE 0. 549 0.00157 1.033 - 1.34337 P »0.10 
HI4(1) HIS(1) 
H1(14) H3(3) Hh (4) 
HH 0.634 0.00186 1.221 0. 38263 P »0.10 
H12(2) HI3(1) a 
H1(10) H2(2) HI2(4) 
JMS 0.618 0.00907 5.956 — 0.64641 P »0.10 
H8(1) 
MDJ H1(12) 0 0 0 - - 
QQHE HI(11) H2(4) 0.419 0.01339 8. 800 1.49070 P »0.10 
SY H1 (11) H2(5) HI7(1) 0.522 0.01435 9.426 1. 78892 0.10 >P >0.05 
SH H1(11)H16(1) 0.167 0.00076 0.500 — 1.62929 0.10 >P 20.05 
RI 
名 种 群 0.456 0.00676 4.440 —0.97226 P>0.10 


Total population 


为 了 进一步 分 析 这 10 个 地 理 种 群 的 遗传 多 样 





性 ,运用 Kimura 2-Parameter 对 已 获得 的 153 条 COI 
基因 序列 进行 了 不 同 单 倍 型 之 间 的 遗传 距离 的 分 析 
( 表 4) ,不 同 单 倍 型 间 遗 传 距 离 为 0. 002 ~ 0.038 之 
间 ,平均 遗传 距离 为 0.012 ,其 中 取 71 H5 单 倍 型 的 
遗传 距离 都 大 于 0. 031 , 均 达 到 了 较 高 的 水 平 。 
2.3 大豆 食心虫 不 同 地 理 种 群 间 的 遗传 多 样 性 、 基 
因 流 及 遗传 分 化 分 析 

利用 固定 系数 Fst 和 基因 流 Nm 对 种 群 间 的 基 
因 差 异 程度 及 交流 程度 (Hudson et al., 1992) 进行 
分 析 。 根 据 DnaSP 5. 0 计算 的 结果 表明 ,总 种 群 的 
Fst 为 0. 12545 , Gst 为 0. 06326 , 总 基因 流 Nm 为 























3.49 , 且 各 种 群 间 的 基因 流 均 大 于 1( 表 5) ,说 明 群 
体 间 基因 交流 的 水 平 较 高 ,通过 某 种 渠道 群体 间 可 
以 发 生 基 因 交 流 , 种 群 间 还 未 产生 明显 的 分 化 。 各 
种 群 间 的 核 背 酸 差 异 数 (Kaxy) 在 0. 1539 ~ 8. 7020 
之 间 , 平 均值 为 4. 3276; 1% E RR BC BE ( Dxy ) 在 
0. 0002 ~0.0132 范围 内 ,平均 值 为 0. 0066 ; 遗传 分 
化 系数 (Cs) 在 — 0. 0252 ~ 0. 1280 之 间 , 其 平均 值 
为 0. 0355 ; fp BER] El XE. RS CC Fst) Æ — 0. 0692 ~ 
0.2512 ,平均 值 为 0.1147。 

应 用 Arlequin 3. 5. 1. 2 软件 进行 AMOVA 分 
子 变异 分 析 ( 表 6) 表 明 ,10 个 地 理 种 群 之 间 的 方差 
组 分 为 0.27252, 占 方差 比率 的 12. 13% ;而 种 群 内 
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X4 大 豆 食心虫 COI 基因 不 同 单 倍 型 间 的 遗传 距离 
Table 4 Genetic distance of different haplotypes of COI gene of Leguminivora glycinivorella 
H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H11 H12 H13 H14 H15 H16 H17 
H1 
H2 0.033 
H3 0.002 0.031 
H4 0.003 0.033 0.002 
H5 0.035 0.002 0.033 0.035 
H6 0.003 0.030 0.005 0.006 0.031 
H7 0.003 0.036 0.005 0.006 0.038 0.006 
H8 0.005 0.031 0.006 0.008 0.033 0.008 0.008 
H9 0.005 0.031 0.006 0.008 0.033 0.008 0.008 0.009 
H10 0.008 0.035 0.009 0.011 0.035 0.011 0.011 0.012 0.009 
H11 0.003 0.036 0.005 0.006 0.038 0.006 0.006 0.008 0.008 0.011 
H12 0.005 0.035 0.003 0.005 0.035 0.008 0.008 0.009 0.009 0.009 0.008 
H13 0.003 0.036 0.005 0.003 0.036 0.006 0.006 0.008 0.008 0.008 0.006 0.005 
H14 0.002 0.035 0.003 0.005 0.036 0.005 0.005 0.006 0.006 0.009 0.005 0.006 0.005 
H15 0.002 0.035 0.003 0.005 0.036 0.005 0.005 0.006 0.006 0.009 0.005 0.006 0.005 0.003 
H16 0.005 0.035 0.006 0.008 0.036 0.008 0.008 0.009 0.003 0.012 0.008 0.009 0.008 0.006 0.006 
H17 0.003 0.033 0.005 0.006 0.034 0.006 0.006 0.008 0.005 0.008 0.006 0.008 0.006 0.005 0.005 0.008 
R5 大 豆 食心虫 不 同 地 理 种 群 间 的 基因 流入 
Table 5 Gene flow ( Nm) between pairwise populations of Leguminivora glycinivorella 
DQ DH GZL HEB HH JMS MDJ QQHE SY SH 
DQ 
DH 20.9592 
GZL 1.9343 2.6368 
HEB 8.5416 10.5375 2.3539 
HH 2.1096 3.8898 2.1867 2.5562 
JMS 4.1773 7.3370 80.1452 5.2013 9.2847 
MDJ inf 17.4856 1. 8697 6.5822 1.6853 3.9723 
QQHE 1. 8969 2.5656 inf 2.3153 2.1205 53.2634 1.8332 
SY 1. 5392 2.0419 inf 1. 8364 1.7272 13.6643 1.4904 inf 
SH inf inf 2.0484 9.7660 2.3752 4.5403 inf 2.0126 1.6204 


inf; 无 限 大 Infinite. 


表 6 


豆 食心虫 10 个 地 理 种 群 线粒体 COI 
基因 的 分 子 变异 分 析 
Table6 Analysis of molecular variance (AMOVA) of 
mitochondrial COI gene in 10 geographic 


populations of Leguminivora glycinivorella 









































变异 来 源 平方 和 方差 组 分 ”变异 百分率 
H BE : 
Source df Sum of Variance Percentage of 
of variation j sequences components variation 
种 群 间 55.043  0.27252Va 12.13 
Among populations 
涯 内 
sig 种 群 内 : 143 282.408  1.97488Vb 87.87 
Within a population 
P4 zi E 
PORA 152 337.451 2.24740 


Total variance 


部 的 方差 组 分 为 1. 97488 , 占 方差 比率 的 87. 87% ， 
种 群 内 部 的 遗传 变异 远大 于 种 群 间 的 遗传 变异 ,说 
明 中 国 东北 地 区 大 豆 食心虫 的 种 群 遗 传 分 化 主要 来 


自 种 群 内 部 ,而 不 是 种 群 

















间 。 


Kimura 2-Paramter 模型 分 析 大 豆 食 心 虫 不 同 地 





理 种 群 遗 传 距离 表明 ( 表 7) ,不 同 地 理 种 群 遗 传 距 
离 在 0 ~ 0. 2909 之 间 。 而 遗传 距离 和 地 理 距 离 的 
Mantel 相关 性 分 析 表 明 ( 图 1) ,各 种 群 之 间 的 遗传 
距离 和 地 理 距 离 之 间 未 呈现 出 显著 相关 性 (r = 
0.210, P 20.913 >0. 05) ,说 明 大 豆 食心虫 不 同 地 
理 种 群 间 的 遗传 距离 不 是 由 地 理 距 离 远近 决定 的 。 

2.4 大豆 食心虫 不 同 地 理 种 群 COI 基因 单 倍 型 系 
统 发 育 分 析 
按照 种 内 遗传 差异 以 选择 亲缘 关系 较 近 的 种 为 
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R7 大豆 食心虫 10 个 地 理 种 群 间 的 地 理 距 离 (kmy) 的 自然 对 数值 (上 三 角 ) 与 遗传 距离 (下 三 角 ) 
Table7 The natural logarithm of geographical distance (km) (above the diagonal) and pairwise genetic distance 
( below the diagonal) among 10 geographic populations of Leguminivora glycinivorella 
DQ DH GZL HEB HH JMS MDJ QQHE SY SH 
DQ 7.76 8.41 7.18 8.83 8.03 9.02 9.32 9.23 6.81 
DH 0. 0239 7.23 7.06 9.35 8.64 8.65 27 9. 12 7.85 
GZL 0.2585 0.1896 8.15 9.63 8.76 8.22 8.45 8.54 8.57 
HEB 0.0585 0.0474 0.2124 9.05 9.18 8.65 8.79 7.80 6.71 
HH 0.2370 0.1285 0. 2287 0.1956 8.32 8.03 8.22 8.98 8.73 
JMS 0.1197 0.0681 0. 0063 0.0961 0.0539 9.40 8.73 9.98 8.16 
MDJ 0.0000 0.0286 0.2674 0.0760 0.2967 0.1259 9.12 9.62 8.32 
QQHE 0. 2636 0.1949 0. 0000 0.2160 0.2358 0.0094 0.2728 9.17 7.95 
SY 0.3249 0.2449 0. 0000 0.2723 0.2895 0.0366 0.3355 0.0000 9.23 
SH 0.0000 0. 0000 0.2441 0.0512 0.2105 0.1101 0. 0000 0.2484 0.3086 
c 
33 
Rd 
3i 
0 0.1 0.2 0.3 0.4 
遗传 距离 Genetic distance 
图 1 大 豆 食心虫 10 个 地 理 种群 的 遗传 距离 
与 地 理 距 离 的 Mantel 检验 
Fig. 1 Mantel test between genetic distance and geographic 


distance among 10 populations of Leguminivora glycinivorella 








外 群 的 原则 ,运用 Kimura. 2-Parameter 模型 邻接 法 
( Neighbor-joining, NJ) , 以 同 科 胡 桃 小 卷 峨 Cydia 
caryana COI 基因 序列 (GenBank 登录 号 :GU088344. 
1) 作为 外 群 ,与 17 个 大 豆 食 心 虫 COI 基因 单 倍 型 
进行 构 树 。 结 果 显 示 ,本 研究 获得 的 17 个 单 倍 型 处 
于 平行 地 位 , 且 都 与 外 群 分 开 ,形成 了 2 个 较为 明显 
的 支 系 , 单 倍 型 HS 和 H2 聚 为 一 支 ,其 余 单 倍 型 聚 
为 一 支 。 尽 管 有 距离 较 近 的 地 理 种 群 的 单 倍 型 聚 为 
一 文 (如 HEB 的 HL 和 DQ B5 H7) , [H8 PR FS E 
远 的 地 理 种 群 的 单 倍 型 聚 为 一 文 ( 如 GZL 的 H5 和 
JMS 的 H2) , 单 倍 型 之 间 未 显示 出 明显 的 地 理 谱系 ， 
各 分 支 置信 度 也 较 低 。 总 体 来 看 ,系统 发 育 树 并 未 
反映 出 与 地 理 位 置 相关 的 信息 ,未 形成 明显 的 系统 
地 理 结构 (图 2)。 

单 倍 型 网 络 中 介 图 能 进一步 说 明 各 单 倍 型 之 间 
的 关系 ,利用 Network 4. 6 程序 构建 单 倍 型 网 络 中 



























































2 邻接 法 构建 大 豆 食心虫 不 同 地 理 
COI 基因 单 倍 型 间 系 统 树 
Fig. 2  Neighbor-joining phylogenetic tree of haplotypes of 





种 群 








mitochondrial COI gene in different geographic populations 


of Leguminivora glycinivorella based on nucleotide sequence 
以 胡桃 小 卷 蛾 为 外 群 。Cydia caryana was used as the outgroup. 

















介 图 (图 3) ,根据 结果 可 知 , 单 倍 型 Hl 与 其 他 单 倍 
型 关系 紧密 ,为 10 个 地 理 种 群 的 共享 单 倍 型 , 单 倍 
型 H2 fü H5 与 其 他 单 倍 型 有 明显 分 支 ,与 发 育 树 结 
果 吻 合 (图 2)。 单 倍 型 的 分 布 与 其 所 在 的 地 理 种 群 
之 间 不 存在 明显 的 对 应 关系 ,大 多 数 单 倍 型 都 处 于 
较为 混杂 的 分 布 格局 中 。 其 结果 更 加 验证 了 2.3 15 
的 分 析 , 说 明 种 群 间 遗 传 分 化 主要 来 自 种 群 内 部 ,种 
群 间 遗 传 分 化 相对 较 小 。 
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本 研究 利用 mtDNA. COI 基因 片段 中 进化 速率 较 
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K3 ”大豆 食 心 虫 线粒体 COI 基因 各 单 倍 型 的 中 介 网 络 图 





Fig. 3 Median-joining haplotype network of Leguminivora glycinivorella based on COI gene of mtDNA 
































圆 面积 代表 单 倍 型 出 现 频率 ,彩色 扇形 面积 代表 各 样品 种 群 在 同一 单 倍 型 中 所 占 的 比例 。Circle areas are proportional to haplotype frequencies, 


while colored portions represent the proportions of the same haplotype that occurs in each sampling region. 


快 的 一 段 序列 作为 分 子 标记 (Bensasson et. al., 
2000) ,对 东北 三 省 大 豆 食 心 虫 危 害 较 重 的 10 个 地 
理 种 群 内 和 种 群 间 的 遗传 多 样 性 进行 了 分 析 。 在 
10 个 地 理 种 群 中 共 获 得 COI 基因 序列 153 条 ,分 成 
17 个 单 倍 型 ,其 中 只 有 单 倍 型 HT 为 共享 单 倍 型 ,说 
明 该 单 倍 型 是 能 够 在 种 群 中 稳定 存在 的 优势 单 倍 型 
( 张 立 勋 等 ,2005 ) ,其 可 能 是 大 豆 食心虫 进化 史上 
出 现 的 较为 稳定 、 环 境 适 应 性 较 强 的 一 种 单 倍 型 。 
而 各 种 群 的 独 有 单 倍 型 说 明 种 群 间 存 在 着 一 定 程度 
的 遗传 分 化 。 在 众多 影响 种 群 遗传 多 样 性 的 因素 
中 , 历史 因素 在 种 群 进化 中 发 挥 很 大 的 作用 (Xiao 
et al., 2009; Yamaguchi et al., 2010) 。 本 研究 中 总 
群体 和 各 种 群 的 Tajima' s D 检验 结果 皆 不 显著 , 表 
明 这 10 个 地 理 种 群 的 大 豆 食 心 虫 在 较 近 的 历史 时 
期 内 并 未 出 现 种 群 扩张 和 持续 增长 的 现象 ,种 群 大 
小 保持 相对 稳定 的 状态 。 

本 研究 中 总 种 群 单 倍 型 多 态度 处 于 中 低 等 水 
平 ,说 明 大 豆 食心虫 对 变化 的 外 界 环境 的 适应 能 
相对 较 弱 ,这 可 能 与 其 生物 学 特性 有 关 。 根 据 李 春 
选 等 (2003 ) 应 用 等 位 酶 的 分 析 可 知 ,昆虫 会 因 其 生 
物 学 特性 及 对 栖息 环境 的 依赖 程度 ,体现 出 不 同 的 
遗传 多 样 性 水 平 。 推 测 大 豆 食心虫 因 其 栖息 环境 、 
食物 的 选择 较为 局 限 ,生存 环境 比较 单一 ,造成 了 大 
豆 食 心 虫 种 群 内 相对 较 低 的 遗传 多 样 性 。 而 单 倍 型 





























多 样 度 和 核 苷 酸 多 样 度 在 各 种 群 中 存在 差异 ,可 能 
是 由 于 大 豆 食 心 虫 在 中 国 东 北 的 分 布 范围 较为 广 
泛 , 所 处 的 自然 生态 条 件 不 同 ,在 自然 选择 的 作用 下 
导致 不 同 地 理 种 群 的 遗传 结构 发 生 改 变 , 而 种 群 的 
遗传 结构 与 物种 内 基因 交流 程度 的 关系 是 重要 的 进 
化 生态 学 问题 ( 祖 元 刚 等 ,1999 ) 。 从 本 研究 结果 来 
看 ,黑河 地 区 种 群 遗 传 多 样 性 最 为 丰富 ,其 次 为 佳 木 
斯 种 群 。 有 研究 表明 ,大 豆 食心虫 属 “ 长 日 照 昆 
虫 ” ,在 每 日 光照 小 于 15 h 出 现 沾 育 ,在 其 盛 发 期 
7 -8 月份 ,高 纬度 的 黑河 地 区 光照 时 间 长 于 低 纬度 
的 哈尔滨 等 地 区 ,因此 ,高 纬度 的 黑河 地 区 大 豆 食 心 
虫 发 生 期 早 于 低 纬 度 的 哈尔滨 等 地 区 ( 张 武 等 ， 
2013 ) 。 由 此 可 以 推测 ,不 同 分 布 区 的 大 豆 食心虫 
发 生 规 律 也 可 能 是 导致 其 遗传 多 样 性 发 生变 化 的 原 
因 之 一 。 另 外 ,本 研究 在 各 地 区 采用 随机 取样 的 方 
法 进 采 集 , 关 于 不 同 品 种 寄主 植物 是 否 会 对 大 豆 食 
心 虫 的 生物 学 特性 等 方面 产生 影响 的 研究 目前 尚 没 
有 报道 ,但 国外 有 研究 表明 母 管 葡 马 属 Mesothrips 的 
葡 马 种 类 在 榕 属 Ficus 不 同 品种 寄主 上 其 行为 特性 
和 生活 史 等 也 不 相同 (Tree and Walter, 2009), PF 
以 基于 上 述 研 究 分 析 , 可 能 本 研究 所 选 的 栽培 品 
由 于 连年 重茬 等 寄主 植物 自身 的 因素 会 对 大 豆 食 心 
虫 种 群 间 遗 传 多 样 性 和 基因 交流 具有 一 定 的 影响 。 
基因 流 Nm 和 遗传 分 化 指数 Fst 可 以 揭示 种 群 
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间 可 能 存在 的 基因 渗透 和 影响 遗传 分 化 的 遗传 现象 
( Rousset, 1997)。 本 研究 的 基因 流 Nm 为 3.49 > 
1, 是 各 种 群 内 部 的 Nm 值 均 大 于 1 ,说明 种 群 间 和 种 
群 内 的 基因 交流 频繁 。 遗 传 分 化 指数 Fst ERIS RR 
传 分 化 不 明显 ,这 与 大 豆 食心虫 群体 间 基 因 交 流程 
度 较 高 的 结论 相对 应 。 此 结果 表明 ,种 群 间 的 基因 
交流 程度 并 不 完全 是 由 昆虫 的 迁 飞 能 力 决 定 的 , 迁 
飞 能 力 弱 或 者 是 无 迁 飞 能 力 的 昆虫 其 种 群 间 也 可 以 
通过 某 些 渠道 进行 基因 交流 。 大 豆 食 心 虫 成 虫 的 飞 
翔 能 力 虽 然 弱 ,但 也 不 排除 其 能 随 气 流 扩散 的 可 能 
性 ,目前 国内 外 还 没有 大 豆 食心虫 成 虫 扩散 等 方面 
的 详细 研究 。 但 是 随 着 外 界 环境 的 压力 和 种 群 结 构 
的 变化 ,种群 间 或 种 群 内 是 否 会 发 生 某 种 遗传 分 化 ， 
有 竺 于 进一步 的 研究 。 一 个 物种 对 环境 变化 的 适应 
能 力主 要 取决 于 其 群体 内 部 的 遗传 多 样 性 水 平和 相 
应 的 遗传 结构 ( 祖 元 刚 等 ,1999)。 从 本 研究 结果 来 
看 ,大 豆 食心虫 的 遗传 变异 主要 来 自 种 群 内 部 而 不 
是 种 群 间 ,说 明 大 豆 食心虫 需要 在 种 群 内 通过 遗传 
变异 来 抵抗 外 界 环境 的 变化 ,使 其 种 群 在 自然 选择 
的 作用 下 能 够 稳定 存在 。 

本 研究 所 分 析 的 10 个 大 豆 食 心 忠 地 理 种 群 , 均 
为 大 豆 食心虫 的 主要 危害 区 及 栽培 大 豆 的 主要 种 植 
产 区 ,地 理 距离 存在 一 定 的 跨度 ,但 相关 性 检验 结果 
说 明 遗 传 距离 与 地 理 距 离 之 间 没 有 显著 线性 相关 
性 。 综 合 上 述 结果 可 知 ,大 豆 食心虫 各 种 群 间 的 基 
因 交 流 并 未 受到 地 理 距 离 的 影响 ,地 理 距离 并 不 能 
解释 种 群 间 的 遗传 变异 。 本 研究 工作 在 一 定 程度 上 
为 大 豆 食心虫 的 系统 进化 研究 提供 了 分 子 水 平 上 的 
理论 依据 ,为 进一步 研究 中 国 大 豆 食 心 虫 的 遗传 多 样 
性 及 其 种 群 分 化 程度 提供 了 理论 和 实验 数据 。 本 研 
究 仅 以 线粒体 COL 基因 作为 分 子 标记 进行 检测 , 且 采 
集 地 点 仅 为 中 国 东 北三 省 。 因 此 , 若 要 全 面 了 解 中 国 
大 豆 食 心 虫 的 遗传 变异 和 地 理 种 群 分 化 程度 及 分 化 
格局 ,还 需要 进一步 扩大 采样 范围 和 样本 量 ; 同 时 采 
用 更 多 的 分 子 标 记 进 行 验证 ,才能 更 加 全 面 地 了 解 大 
豆 食心虫 各 分 布地 区 的 种 群 间 关 系 和 遗传 机 制 。 
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